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Open issues in stellar modelling 
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ABSTRACT 

An important goal of helio- and asteroseismology is to improve the modelling of stellar evolution. Here I provide a brief 
discussion of some of the uncertain issues in stellar modelling, of possible relevance to asteroseismic inferences. 
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1. Introduction 

The goal of the workshop was to investigate ways to im- 
prove our understanding of the Sun; this is obviously in- 
timately linked to the general understanding of stellar 
structure and evolution, and indeed there are considerable 
prospects that our growing possibilities of asteroseismic 
sounding of other stars will inform our studies of the so- 
lar interior. 

From the point of view of asteroseismology, the relevant 
aspects of stellar modelling include both the study of stellar 
structure and evolution and the modelling of stellar oscil- 
lations, in particular their frequencies, for a given model. 
The latter aspect provides the diagnostic link between the 
observations and the stellar models; although the adiabatic 
approximation is valid for the oscillations in most of the 
star, departures from adiabaticity, and other uncertainties 
in the modelling of the near-surface layers, give rise to sub- 
stantial systematic errors that must be taken into account 
in the analyses. 

Unlike what is perhaps a common perception, we 
are still far from adequate modelling of stellar interi- 
ors. Here I can only touch on a few issues, mainly in 
connection with the modelling of main-sequence stars 
showing solar-like oscillations. A more detailed discussion 
of these issues, and further references, was provided by 
IChristensen-Dalsgaard fc Houdekl (|2010D . For an extensive 
presentation of stellar oscillatio ns and helio- and a steroseis- 
mic techniques and results, see lAerts et al.l (|2009D . 

2. Numerical issues 

A prerequisite for meaningful asteroseismic diagnostics of 
the physics of stellar interiors is that stellar modelling 
presents a faithful representation of the physical assump- 
tion. Thus the models must be numerically sufficiently ac- 
curate. Perhaps the best test of this is to compare indepen- 
dently computed models with the same physical assump- 
tions. A major effort towards such comparisons was carried 
out in the ESTA project initiated as part of the prepara- 
tions for analysis of CoRoT data (|Monteiroll2008( ). In gen- 
eral, the agreement was reasonable between the results of 
the evolution codes included in the comparison, although 
not obviously adequate for the asteroseismic analysis; also, 
it would be very valuable to extend the comparison to other 



codes commonly used for general stellar modelling. The re- 
sults of adiabatic oscillation calculations for a given model 
agreed quite well, but the analysis highlighted the impor- 
tance of adequate numerical resolution in the evolution and 
oscillation calculations, and of consistency in the descrip- 
tion of the stellar physics. 

An important, and probably often inadequately treated, 
aspect of the computations is the accurate specification, 
implementation and documentation of the unavoidable ap- 
proximations in stellar modelling. Without appropriate 
care in this area it is difficult or impossible to use aster- 
oseismic inferences to test the validity of these approxima- 
tions. 



3. Stellar parameters 

Efficient utilization of the asteroseismic data requires the 
best possible information about other properties of the star. 
In the solar case, the mass, radius, luminosity and age are 
determined quite accurately (or at least precisely) from in- 
dependent observations. Solar composition, characterized 
by the ratios of abundances of elements heavier than he- 
lium to the abundance of hydrogen, can be determined 
from spectroscopic observations. Recently, however, there 
has been a substantial revision in some of these abundances, 
leading to conflicting comparisons between t he resulting so- 
lar m odels and helioseismic inferences (e.g.. lAsplund et al.l 
l2009l and references therein). 

Parameters of other stars are in general known far less 
well. Masses can be obtained in the rare cases where the 
star is a member of a well-observed binary system. The 
effective temperature and surface gravity can be obtained 
from spectroscopic observations, but subject to the possible 
limitations of modelling of stellar atmospheres and hence 
with substantial (and probably often underestimated) un- 
certainties. The stellar luminosity requires knowledge, from 
parallax observations, of the distance and hence is cur- 
rently restricted to relatively nearby stars; also, to the ob- 
served stellar magnitude must be applied a bolometric cor- 
rection which again depends on atmosphere models. Also, 
in a few cases the stellar radius can be determined from 
interferometry, again assuming that the distance is known 
and with some sensitivity to atmospheric structure through 
limb darkening. Stellar composition is obtained from spec- 
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troscopy; in the case of stars similar to the Sun this is most 
often done differentiaUy, relative to the solar spectrum, and 
hence the abundances are directly affected by the uncer- 
tainty in the solar composition. 



4. Microphysics 

For the purpose of stellar modelling the equation of state 
is probably in general sufficiently well known from re- 
cent tabulations of sophisticated equations of state. On 
the other hand, helioseismology clearly demonstrates that 
these are not yet co rrect , a t the level of the observational 
precision (e.g., [Basu et all ll999D . This would also affect 
the asteroseismic determination of the helium abundance 
from the signatures in th e frequencies of helium ionization 
(|Houdek fc Gough|[2007h . 

The computation of stellar opacities is considerably 
more uncertain than the equation of state, with direct 
effect on the structure of the radiative part of stellar 
models. Since the heavy-element abundance affects stel- 
lar structure predominantly through the opacity, uncer- 
tainties in the heavy-element abundances and the opac- 
ity are closely linked. Thus an obvious way to correct 
solar models, given the revised abui idances, is to claim 
substantial opacity increases (e. g., iBahcall et al.l l2005l 
IChristensen-Dalsgaard et al] I2009D . although possibly be- 
yond what is physically realistic. An independent indica- 
tion of a need for opacity increases, although at somewhat 
lower temperatures than relevant in the Sun, comes from 
the lack of pred icted instability of some observed m odes in 
P Cephei stars ()Dziembowski fc Pamvatnvkhll2008[ ). 

Although there remain substantial uncertainties in im- 
portant nuclear parameters the effect on stellar modelling 
is in general relatively modest, since even large changes in 
the parameters can be compensated by modest changes in 
the temperature, owing to the high temperature sensitiv- 
ity of the reactions. An important exception concerns the 
balance between contributions to the PP chains and the 
CNO cycle in hydrogen burning, which has a substantial 
effect on the presence and extent of convective cores. This 
includes the relatively recent large reducti on in the rate of 
proton capture by ^'^N (jAngulo et al.|[2005l ). An interesting, 
and so far not res olved, issue concerns electron s creening of 
nuclear reactions (|Shavivll200i iMao et al.ll2009t) . 

There is no doubt that diffusion and settling take place 
in those parts of a star where there is no macroscopic mo- 
tion. These processes, in various approximations, are now 
universally included in 'standard' solar modelling, leading 
to a increase of a few per cent in the surface hydrogen 
abundance during evolution to the present solar age, and a 
decrease of around 10 per cent in the heavy-element abun- 
dances. In somewhat more massive stars with thinner outer 
convection zones the settling rate at the base of the con- 
vection zone is much higher, leading to an almost complete 
elimination at the stellar surface of helium and heavier el- 
ements, o n a timescale short c ompared with the evolution 
timescale (jVauclair et al.l[r974l ). To account for the 'normal' 
abundances observed in most such stars one must therefore 
invoke mixing processes or possibly mass loss to compensate 
for the settling. A possible explanation is mixing caused by 
rotationally induced meridional circulation (sec below). 

An additional complication, particularly in stars a lit- 
tle more massive than the Sun, is the selective effects of 



radiation pressure on different elements, leading to gravi- 
tational levitation counteracting settling and strong local 
variations in the heavy-element composition. To be taken 
properly into account, this requires opacity calculations for 
the local composition, depending on loca tion and time, as 
the star evolves (e.g., [Richer et al.ll200dl ). This has so far 
only been consistently implemented in very few evolution 
calculations. 

5. Properties of convective cores 

Convective cores play an important role in the main- 
sequence evolution of stars of masses just slightly higher 
than the Sun and above. This is caused by the increasing 
dominance in hydrogen burning of the much more temper- 
ature sensitive CNO cycle over the PP chains. For stars 
of masses less than around 2 solar masses this involves 
a phase where the mass of the convective core increases, 
owing to the gradual conversion, on a timescale compara- 
ble to the evolution timescale, of ^^O to ^''N; if diffusion 
is neglected the growth of the core leads to a discontinu- 
ity in the hydrogen abundance at the edge of the core, 
and hence in the density and sound speed. Asteroseismic 
diagnostics of this discontinuity, and other aspects of 
convectiv e cores, may be possible with sufficiently accu- 
rate data (iPopielski &: D zicmb owskil2005llMazumdar et al] 
l2006LICunha fc Metcalfell2007l) . 

The uncertainties in the microphysics, particu- 
larly as it affects the importance of the CNO cy- 
cle, influence the size of the con vective core (see 
IChristensen-Dalsgaard fc Houdekl [20101 for details). Thus 
the reduction in the ^'^N reaction rate shifts the onset of 
convective cores higher in stellar mass by about 0.06 Mq. 
An interesting case is the effect of the revision of so- 
lar abundances which, as discus sed above, is reflected in 
the assumed stellar abundances. IVandenBerg et all (|2007l ) 
showed that this led to a significant change in the isochrones 
computed for the open cluster M67; with the old composi- 
tion models near the end of the central hydrogen burning 
had a convective core, as also suggested by the observed 
colour-magnitude diagram, while models with the revised 
composition lacked the convective core. 

A probably more important uncertainty concerns the 
extent of convective overshoot. There is little doubt that 
motion continues beyond the convectively unstable region, 
but the extent of that motion, and its effects on stellar 
structure, are highly uncertain. Presumably the motion is 
sufficiently vigorous to cause homogenization of the compo- 
sition, but it is less clear whether it leads to full mixing of 
entropy and hence an adiabatic stratification. The extent 
is typically parameterized as a fraction aov of the pressure 
scale height at the edge of the core, with a correction for 
very small cores, but no a priori estimate of aov is available. 
Analyses of open clusters and binary stars lead to values 
of aov of typically around 0.1 - 0.2, a value c onfirmed by 
asteroseismic analyses of /3 Cephei stars fe.g.. lAerts et al.l 
2003). 

An additional complication in models with growing con- 
vective cores is the presence of what has been called semi- 
convection. Convective instability is typically defined in 
terms of the temperature gradient V — dlnT/dlnp, where 
T is temperature and p is pressure; convective instability 
sets in where the value Vrad of V required to transport 
energy by radiation exceeds the adiabatic value Vad- Since 
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Fig. 1. Properties near the edge of the convective core in 
models of a 1.3 M0 star, as a function of distance to the 
centre in units of the stellar radius; solid lines show a model 
without overshoot, dashed lines a model with modest over- 
shoot, Q!ov = 0.05. The upper panel shows the hydrogen 
abundance by mass, illustrating the discontinuity at the 
edge of the mixed core. The lower panel shows the depar- 
ture of the radiative temperature gradient Vrad from the 
corresponding adiabatic gradient. 

Vrad is proportional to the opacity which increases with 
increasing hydrogen abundance X, Vrad jumps to a larger 
value at the edge of a growing convective core where X is 
discontinuous. This is illustrated in Fig. 1, in a model calcu- 
lation where the extent of the convective core is defined with 
marginal instability evaluated for the composition outside 
the core, leading to a small convectively stable region in the 
outer parts of the core. Treating the mixed region in this 
manner is obviously a computational artifice, rather than 
a physically justified model. Modelling of this region, with 
added complications when the settling of heavy elements is 
taken into account, tends to lead to irregular fluctuat ions in 
the m ass contained in the convective core (Lcbrcto n et al.l 
I2008f) . As illustrated, even a modest amount of overshoot 
shifts the composition discontinuity sufficiently far into the 
stable region that the problem is eliminated. 

6. Near-surface problems 

The treatment of convective envelopes also involves sub- 
stantial uncertainties. Overshoot below the convective en- 
velope has a relatively modest effect on stellar evolution 
although it can affect the properties of the red bump on 



the red-giant branch. A more serious concern are the prop- 
erties of the near-surface layer where the density is low 
and consequently a substantial superadiabatic gradient is 
required to transport the energy. Together with the struc- 
ture of the stellar atmosphere this determines the spe- 
cific entropy in the predominantly adiabatic bulk of the 
convective envelope and hence its structure, including its 
depth. In the solar case the t reatment of this layer, e.g., 
usmg the lBohm-Vitensd (|1958D mixing-length formulation, 
is calibrated to obtain the correct radius; this calibration is 
typically, with little justification, used in modelling other 
stars. Hydrodynamica l simulations of near-surface convec- 
tion (|Nordlund et al.l 12009) provide a reasonably realistic 
modelling of these layers; unlike other parts of the star the 
relevant dynamical and thermal timescales are sufficiently 
similar that the relevant effects can be taken into account, 
although obviously still with an approximate treatment of 
scales smaller than the numerical resolution. The results of 
the simulations ca n then be used to calibrate the simpler 
formulations fe.g.. lTramt)edachll2007f) : this offers a promis- 
ing procedure for more realistic modelling of this part of 
the star, although it has so far not seen much use. 

Uncertainties in the modelling of the near-surface layers 
have a substantial effect on the oscillation frequencies and 
hence on their use as asteroseismic diagnostics. In addition 
to the structure of the superadiabatic layer, these uncer- 
tainties include the dynamical effects, usually ignored, of 
convection on stellar structure in the form of 'turbulent 
pressure', nonadiabatic effects on the oscillations, and the 
coupling between convection and pulsations, in terms of 
the perturbation to the convective flux and the turbulent 
pressure, as well as the stochastic excitation of the modes, 
for solar-like oscillations. These effects dominate the dif- 
ference between the observed and modelled frequencies of 
solar oscillations; they can be suppressed, however, in he- 
lioseismic analyses because of the broad range of degrees of 
the observed modes. In the stellar case this is not possible, 
i n gen eral. It was pointed out by [Roxburgh fc Vorontsovl 
(j2003l ) that combinations of frequency separations can be 
constructed which are insensitive to the superficial layers 
and retain their sensitivity to the properties of the core (see 
also lOti Floranes et al.l[2005l ). For more general use of the 
frequencies, including calibrations of the overall properties 
of the star, one can attempt to estimate the near-surface 
effects on the frequencies, by assuming a f unctional form 
simil ar to the known effect in the solar case (jKieldsen et al.l 
I2OO8D . This, however, remains a serious issue in asteroseis- 
mic analyses. 

7. Rotation 

There is no doubt that most, or indeed all, stars rotate, yet 
rotation is usually ignored in modelling of stellar evolution. 
A detailed discussio n of the effects of rotation on stars was 
recently provided bv iMaedeil ()2009D . 

The dynamical effects of rotation on stellar structure 
are relatively straightforward to incorporate, at least as 
long as they can be treated as perturbations around a non- 
rotating, spherically symmetric structure. In a slowly ro- 
tating star such as the Sun these effects are very small. 
However, many stars rotate so rapidly that the pertur- 
bative approach is inadequate; here two-d imensional mod- 
elling of stellar structu re is required fe.g.. lRoxburgh|[2bo4 
iMacGregor et al.ll2007D . Far greater comphcations are as- 
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sociated with the effects on stellar evolution, including the 
evolution of the internal rotation rate. A naive local ap- 
plication of the conservation of angular momentum would 
predict that the angular velocity of the central parts of stars 
will increase with age as these regions contract, while rota- 
tion in the outer parts would be expected to slow down as 
they expand. This is certainly too simple. As already noted 
by von Zeibel and Eddington, rotation causes a thermal im- 
balance which leads to circulation and hence redistribution 
of angular momentum and mixing of the stellar composi- 
tion. Indeed, it is likely that in many stars this mixing coun- 
teracts the rapid settling discussed above. To these pro- 
cesses must also be added mass loss, possibly magnetically 
linked to the stellar convective envelope, which removes an- 
gular momentum from the star. It seems likely that most 
stars start their life with rapid rotation; stars with masses 
up to somewhat higher than the Sun apparently lose angu- 
lar momentum to a magnetized stellar w ind, leading t o a 
strong decrease in rotation with age (e.g.. lBarnesll2003[ ). 

A treatment of these processe s was proposed by IZahnI 
(|l992f) and further developed by iMaeder fc ZahnI (|l998t) . 
This assumes an angular velocity that depends only on 
the distance to the centre of the star, as a result of strong 
horizontal turbulence. Mixing of composition is a diffusive 
process while transport of angular momentum in addition 
includes advective terms. This formulation has seen fairly 
extensive use and has had some success in accounting for 
the observed composition of massive stars. 

A serious problem is to account for the helioseismi- 
cally inferred solar internal rotation rate (e.g. . Illow3l2009t ) : 
in particular, the Zahn model is unable to explain the 
present slow rotation of the radiative interior. This re- 
quires additional mechanisms transporting angular mo- 
mentum from the interior to the convection zone. It has 
been proposed that this coupling could be mediated by 
gravity waves excit ed at the base of the convection zone 
(iMathis et al. 20081): alternatively, it may be of magnetic 
nature ( Garaud fc Guervillvl 120091 ). It is obvious that as- 
teroseismic information about the internal rotation of other 
stars, although unavoidably quite limited in the foreseeable 
future, can be extremely valuable in distinguishing between 
these mechanisms. 



8. Concluding remarks 

It is evident that there are many serious open issues in stel- 
lar modelling. An important task is the evaluation of the 
asteroseismic signatures of these effects, including the de- 
sign of diagnostics that may best investigate them and a 
determination of the resulting requirements on the obser- 
vations. The asteroseismic observations that are currently 
been obtained by the CoRoT and Kepler space missions 
certainly provide excellent prospects for addressing these is- 
sues, although the experience from CoRoT has shown that 
the analysis of the data also involves serious challenges. 

There is clearly a need for very substantial development 
of the techniques of stellar modelling. This can be inspired, 
but certainly not replaced, by further detailed hydrodynam- 
ical simulations of specific aspects of stellar interior dynam- 
ics. Except for the near-surface layers a serious constraint 
is the huge mismatch between the relevant dynamical and 
thermal timescales, implying that the simulations cannot 
be run under realistic stellar conditions. A great deal of 



physical insight will therefore be required to extrapolate 
the results of the simulations to those conditions. 

These efforts will require intensive collaborations be- 
tween data analysis, data interpretation, theory and mod- 
elling. A workshop such as the present is an ideal venue for 
furthering such collaboration and create new ideas. 
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